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I present results for two-loop soft anomalous dimensions, which are derived from dimensionally regularized 
diagrams with eikonal quark lines and control soft-gluon emission in hard-scattering processes. Detailed results 
for the UV poles of the eikonal integrals are shown for massive quarks, and the massless limit is also taken. The 
construction of soft anomalous dimensions at two-loops allows soft-gluon resummations at NNLL accuracy. 
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1. Introduction 

Soft-gluon resummation formalisms for top quark 
production and other hard scattering cross sections 
beyond leading logarithms are formulated in terms of 
exponentials of soft anomalous dimensions that con- 
trol noncollinear soft-gluon emission (see e.g. Gil- 
The calculations of these anomalous dimensions [2, 3] 
are performed using the eikonal approximation, which 
is valid for descibing the emission of soft gluons from 
partons in the hard scattering. The approximation 
leads to a simplified form of the Feynman rules by 
removing the Dirac matrices from the calculation. 
When the gluon momentum, k, goes to zero, the 
Feynman rule for the quark-gluon vertex reduces to 
gsTpV^ jv ■ k with v a dimensionless velocity vector. 

Here we calculate diagrams and derive the soft 
anomalous dimension for top quark pair production 
via e + e~ — ► ti through two loops 3. We also discuss 
the massless limit and extensions to other processes. 

Writing the soft anomalous dimension Ts as a series 
in a s 



7T 



a s Vr< 2 ) 



(1) 



we calculate the one- and two-loop expressions for Ts ■ 
We use the Feynman gauge and we employ dimen- 
sional regularization with n = 4 — e dimensions. The 
soft anomalous dimension can be determined from the 
coefficients of the ultraviolet (UV) poles of the eikonal 
diagrams. 

We do not include the color factors in the individ- 
ual diagrams below, but take them into account when 
assembling all the pieces together in the last section. 
The Appendix lists some of the integrals calculated 
for the evaluation of the diagrams. 



2. One-loop diagrams 

In this section we calculate the one-loop diagrams 
in Fig. 1. We define — v/l — 4m 2 /s with m the 
heavy quark mass and s the squared cm. energy. We 
label the two heavy quark lines by i and j. Using the 





(a) 



(1.) 



Figure 1: One- loop diagrams with heavy-quark eikonal 
lines. 



relation v — p W2/s, with p the momentum, we have 
v 2 = v 2 = (1 - /3 2 )/2 and m ■ Vj = (1 + /3 2 )/2. 

We begin with the integral I\ a for the diagram in 
Fig. 1(a) given by 



ha = g s 



d n k (-i) 5p „ < (-t#) 



(27T)™ k 2 



k (-Vj ■ k) 



(2) 



Using Feynman parameterization and after several 
manipulations and separation of UV and infrared (IR) 
singularities, we find through Eq. (|69|) the UV poles 



ha = 



a s (l + /5 2 )l (1-0 



20 



In 



1+0 



(3) 



Next we calculate the one-loop self-energy diagrams 
in Fig. 1(b) given by 



hb = 9t 



d n k (-i)g^ 



(27t)" k 2 



[k 1 -k) Vi- k' 



(4) 



Here we have introduced the regulator k', i.e. the 
external quark momentum is pt + k', so as to use the 
eikonal rules. We then expand the above expression 
around vi-k' = at constant e. The expansion gives an 
irrelevant 1/vi ■ k' term, a constant term, and linear 
and higher-order terms in «j • k' which vanish when 
setting k' = 0. The remaining term is then 



'16 



i9 2 s v 2 



d n k 



(2tt)™ J k 2 (Vi ■ k) 2 ' 
We isolate the UV poles and using Eq. (|7D|) we find 



(5) 



Ct s 1 
lib = • 

7r e 



(G) 
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Figure 2: Two-loop vertex diagrams with heavy-quark 
eikonal lines. 



The one-loop soft anomalous dimension is then read 
off the coefficient of the UV poles of the one-loop di- 
agrams through 



Cf[Iici + hb\ 



which gives 



W) h/ 1 



2/3 



■0 



1 + 



(7) 



(8) 



with C F = (N 2 - 1)/(2JV C ) the color factor. 

When expressed in terms of the cusp angle 7 = 

cosh" 1 fvi ■ Vj/^Jvfvf] = ln[(l + - /?)] this be- 

comes = Cp (700th 7 — 1) and is also known as a 
cusp anomalous dimension [4|, |5[ . 



3. Two-loop vertex diagrams 

The two-loop diagrams can be separated into two 
groups, shown in Fig. 2 and 3. The calculations are 
challenging due to the presence of a heavy quark mass 
and they involve multiple complicated integrals and 
delicate separations of infrared and ultraviolet poles, 
which by construction of the soft function are oppo- 
sites of each other Q . The analytical structure of the 
results involves logarithms and poly logarithms. 

In this section we calculate the two-loop vertex dia- 
grams of Fig. 2. We also have to include the one- loop 
counterterms to these diagrams. We will denote by I 2 
the diagrams in Fig. 2 and by 1% their counterterms. 
The total contribution will then be I 2 = I 2 + 1% ' '• 

We find it convenient to define 



Mp = (41n2 + ln7r — 7b — i7r) In 



1-/3 
1 + 



+ iln 2 (l + /3)-iln 2 (l-/3) 



Li 2 



1+0 



Li 2 



1-0 



(9) 



3.1. 1 2a and I 



2b 



In this section we calculate the integrals I 2a and 
I 2 b for the two-loop diagrams in Figs. 2(a) and 2(b). 
Note that there is no counterterm for diagram 2(b), 
so I 2b = V 2b . 

Diagram 2(b) is given by 



r 4 

hb = 9 S 



4 f d n ki d n k 2 {-i)g^ v (-i)g p 



(2tt)" (2tt)™ k{ 



k 2 



{-<) (-<) 



Vi ■ ki Vi ■ (fci + k2) —Vj ■ (fci + k 2 ) —Vj ■ k 2 



(10) 



Now, we begin with the k 2 integral and use Feyn- 
man parameterization. Performing first the k 2 inte- 
gration and then an integration over one of the pa- 
rameters, we have 



i 2b = - A 2- 4+e 7T- 2 +^ r ( 1 - 1 ) r(i + e) 



x(l + /3 2 ) 2 / dz / dy(l-yY 
Jo Jo 



2(3 2 (\-y) 2 z 2 -20 2 (\-y)z 



(1-P 2 ) 



k\ Vi ■ fci [((Vi - Vj)z + vj) ■ fci 
Next we proceed with the fci integral 



-1+4 



(11) 



k\ Vi ■ fci [((Vi — Vj)z + Vj) ■ k 



x2 2+3e 7T 2 -* 



r(i + e) Jo 



il+e 



dX!Xi 1+2e (l- Xi) 



-i(-l)" 



-l-2e 



x / dx 2 (l - x 2 y {[x 2 Vi 

JO 



+ (1 - X 2 ) ((Vi - Vj)z + Vj)] 2 j 



-1-34 



(12) 



The integral over x\ contains both UV and IR singu- 
larities. We isolate the UV singularities via 



dx 1 x^ 1+Ae (l-x) 



-l-2e 



dx\ x 1 



l+2e 



+ / d Xl xr 1+2e [(1 - x)- 1 - 2 ' - 1] = 1 + IR. (f 3) 
Jo 2e 

We set e = in the integral over x 2 since it is UV 
finite. We thus find the UV poles of the fci integral 



cTfci 



fcf Vi ■ fci [((Vi - Vj)z + Vj) ■ fci 



= -i 2tt 2 



,1 



1 



0(1 -z) 



[tajah _1 (/3(l - 2z)) - tantT 1 (-/?)] .(14) 
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Using this result in the expression for I 2 b, Eq. (fTTj) , 
and performing the y and z integrals, we find 



'2fc 



a 2 (l + (3 2 ) 2 l 



8/3 2 



hi 



+ Li 3 



1-/3 
1 + /3 

+ 



1 

Li 2 



_/3 
•/3 
(1-/3) 2 
(1 + /?) 2 



(15) 



Diagrams 2(a) and 2(b) are related through the 
equation I 2a + lib = \l\ a (where we need to keep 
UV poles and constant terms in Ii a ) from which we 
can calculate I 2a ■ The one-loop counterterm to I 2a is 



jC.t. 

1 2a — 



* 2 (1+/? 2 ) 2 

7T 2 8/3 2 



2 , 2 



1-/3 
■/J 



>w) b (ir| 



(16) 



Next we calculate the integral for the gluon-loop 
diagram given by 



d n k d n l 



V; 



r =V 

2c 9 i 2 y. j ^27r) n (2n) n Vi ■ k (—Vj ■ k) 



(-j)g^> {-i)9 PP ' M).<w (.-i)g vv > 

k 2 I 2 (k-l) 2 k 2 

g pp (k + l) a + g pa (k - 2iy' + g ap (-2k + l) p 

g p ' u \l + ky' + g"' a ' {-2k + l) p ' 



+ g ap (k-2l) v 



Using Eq. ([71]) we find for the UV poles of I' 2 



(22) 



2cgl' 



_ a 2 19 (l + /3 2 ) 
2cgl 7T 2 96 (3 

1 

e 



1 + /3 



57 V 1 + P 



Mr, 



(23) 



Then we find the simple relation 
a 2 s (l+P 2 ) 2 (-1) 



3.2. / 



/26 = 



8/3 2 



hr 



1-/3 
1 + /3 



(17) 



2 c 



In this section we calculate the integrals for the 
three diagrams represented by Fig. 2(c). The blob 
represents a quark, gluon, or ghost loop. Note that 
an additional diagram with a four-gluon vertex van- 
ishes. 

The quark-loop diagram is given by 



Acq = {- l ) n S9 4 s 



d n k d n l 



(-<) 



icq v / (27r) n (27r) n Vi ■ k (-Vj ■ k) 

(-i)9pv (-i)g P a 



k 2 



k 2 



-Tr 



I 2 



(l-kf 



.(18) 



After several steps, and using Eq. (f7Tj) . the final result 
for the UV poles of I' 2cq is 



Acq — — f n f 



a 2 (1 + /3 2 ) 



6/3 



1 



2 In 



5 m 

6 



1-/3 
1 + /3 
1-/3 
1 + /3 



(19) 



with the number of light quark flavors. The one- 
loop counterterm is 



TC.t. 

1 2cq 



v 2 



n; (1 + /3 2 ) 



6/3 



In 



1-/3 
1 + /3 



-Mr- 



Then the sum / 2 c 9 = Acq + Acq & ves 



hcq — 



v 2 



a: (1+/3 2 ) 



6/3 



5 

67 



ln 



1-/3 
1 + /3 



(20) 
(21) 



The one-loop counterterm is 



t . _a 2 19 (l + /3 2 ) f 2 fl-/3 

^ " 7T 2 96 1 6 2 



1 + /3 



(24) 



Then 



Hcgl 



a 2 19 (l + /3 2 ) f 1 
^96 /T~ 



57eJ V1 + /3 



(25) 

Last we calculate the integral for the ghost-loop di- 
agram given by 

I^ = {-l)gi 1 - - — — 



L 2cgh 



d n k d n l v p (-v p ) 
(2n) n (2n) n Vi ■ k (-vj ■ k) 



x4r- 



i 2 (i - ky 

Using Eq. (J7TJ), we find 



(I - k) 



a {~i)9pu (-i)g pa 



k 2 



k 2 



(26) 



1 2cgh 



(1 + /3 2 ) 
Tr 2 96/3 



1 + /3 



The one-loop counterterm is 



« 2 (l + /3 2 ) 



2c9h Tr 2 96 /3 
Then 

a 2 (l + /? 2 ) f 1 



1-/3 
1 + /3 



^2cg/i — 



7T 2 96^ 



i} ln (rrl 



(27) 



(28) 



(29) 



The sum of the gluon and ghost loops, Eqs. |25|) and 
(|2"9"|) , denoted by i2 CS , is then 



« 2 5 (1 + /3 2 ) 
Tr 2 24 /? 



31 
We 



hi 



1-/T 
1 + /3, 



(30) 
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3.3. hd and I 



2c 



In this section we calculate the integral I 2d for the 
diagram in Fig. 2(d) given by 



I 2d 



d n k d n l (-i)g^ v {-i)g P o 
(2tt)™ (2tt)« k 2 I 2 

Vi ■ (k — I) Vi ■ k {—Vj ■ k) 



Vi ■ k 



(31) 



as well as the integral I 2e for the diagram in Fig. 2(e) 
given by 



4 



9, 



d n k d n l {-i)g MV {-i)g P a 



(2tt)« (2tt)" k 2 



I 2 {-Vi ■ I) 



(-«?) 



Vi ■ (k — l) Vi ■ k (—Vj ■ k) 
First we note that 



(32) 






(a) 



OH 



(c) 



Figure 3: Two- loop heavy-quark self-energy diagrams 
with eikonal lines. 



(3ali) 



(3alii) 



(3a2) 



Figure 4: Detail of the diagrams of Fig. 3(a) 



a.: 



d n k 



{2Tr) 2nVl Vl "° j J k 2 {v t -k) 2 v r k 
d n l 







(33) 



where the last integral is zero because it is odd in 
I. Therefore I' 2d = —I 2e and similarly for the coun- 
terterms, and then I 2d = —he- Thus we only have 
to calculate I' 2e and its counterterm. Using Eq. ([72")) 
and including the counterterm I 2 *- = (a 2 /7r 2 )[(l + 
/3 2 )/(4/3)][2/e 2 ln((l - /?)/(! + 0)) + M /e] we find 



he = 



7T 2 4/3 

a. l ^ 2 
+ -ln 



1-/3 

■J9 



In 



1 + 



1-/9 



- 2 Ll2 i(TW 



In 



1 + /3 



1 + /3 



In 



4/3 



2 



(34) 



2( 2 + In' 



1-/3 
1 + /3 



1 + 



2/3 



(36) 



4. Two-loop heavy-quark self-energy 
diagrams 

In this section we calculate the diagrams in Fig. 3 
and their counterterms. We will denote by I' 3 the dia- 
grams shown in Fig. 3 and by 7| their counterterms. 
The total contribution is then 1^=1^ + 7|*\ As for 
diagram 1(b) we introduce a regulator k' in the quark 
momentum and then expand around Vi • k' = at 
constant e. We also find it convient to define 



K 



13 



-ln(l-/3 2 ) + 51n2 + ln7r- 7jB -i7r. (37) 



3.4. 7 2/ 

In this section we calculate the three-gluon diagram 
in Fig. 2(f) given by 

w _ „4 / d " fc 2 {-i)9nn' {-i)9w> {-i)g PP 



hf = 9s 



(2?r)" (2tt)™ fcf 

? (-<) (-0 



fcf (fci + k 2 ) 2 



Vi ■ ki (—Vj ■ fci) — Uj • (ki + k 2 ) 



H) 



g^(h-k 2 r +g vp (h+2k 2 y 
+ g p '»'(-2k 1 -k 2 r' 
After many steps we find 



1 / (1 + /3 2 ) , 3 fl-0 

l2f = -e\-^r ln {— 



(35) 



4.1. h a 

Here we calculate the three diagrams represented 
by Fig. 3(a) which are shown in detail in Fig. 4. Note 
that an additional graph involving a three-gluon ver- 
tex with all three gluons attached to the same eikonal 
line vanishes. 

We begin with the diagram labeled (3ali). Using 
Eq. ||7D|) we find 



hali 



1 1 



-K, 



(38) 



with counterterm 7^^ 
Then 



halt — 



(a 2 A 2 )[-2/e 2 - K p /e] 



a 2 (-1) 



(39) 
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Next we calculate I^m- Using Eq. (|75|) we find 



(40) 



with counterterm 7| a * Hi = (a 2 / 'tt 2 )[-1 / 'e 2 - Kp/(2e) 
Then 



1 



1 



ham- l\ 2e 2 + 2e 



(41) 



Adding Eqs. (|39[) and ([4"Tj) and denoting the sum by 
I 3a i, we find 



<K 



3 1 
Ti 2 " + 27 



(42) 



Last we calculate I^ a2l given by 
4 /" eTfci d"fc 2 wf 



3a2 — 9s 



(2tt)™ (2tt)™ (fc' - fci) 



Vi ■ (k' — ki — k 2 ) Vi ■ (k' - k 2 ) k{ 



k\ 



(43) 



Using Eqs. ([74). ijTS]). and ([76]). we find the UV poles 
of the integral 



ha2 - n 2\ ,2 e 



(44) 



with counterterm If^ — ( a s/ 7r2 )[2/e 2 + Kp/e]. Then 



'3a2 



(45) 



4.2. 1 3b 

In this section we calculate the three self-energy di- 
agrams (with quark, gluon, and ghost loops) repre- 
sented by Fig. 3(b). Note that an additional diagram 
with a four-gluon vertex vanishes. 

After several manipulations the quark- loop diagram 
becomes 

gj 1 f d n kd n l 
3bq Uf (2tt) 2 ™ Vi-k' J k 4 Vi-(k- k') I 2 (I - k) 2 
x [2(vi ■ I) 2 - v 2 I 2 - 2v t ■ I v t ■ k + I- k] . (46) 

Expanding around v%- k' = this gives 

r lr „ A f d-k r <ri 

3b? »/ ( 2 TT) 2n J k 4 {Vi ■ k) 2 J I 2 {I- k) 2 

x [-v 2 l-k- 2(vi ■ I) 2 + 2v l -lv l -k] . (47) 

A calculation of the UV poles of the integral using Eq. 
([77)) gives 



hbg - n2 3 



1 1 

72 ~ 7 



(48) 



with counterterm = (a 2 /7r 2 )(n / /3)[2/e 2 + Kp/e] 
Then 



hbq 



a 2 n f 


' 1 


5 


IT 2 3 


.e 2 


6e 



(49) 



Next we calculate the gluon-loop self-energy dia- 
gram in Fig. 3(b) which, after several manipulations, 
becomes 

Jl - 



i 



d n k 



2(2ir) 2n Vi ■ k' J Vi-(k-k') 



d n l 



- (An - 6) 
(An - 6) 



Vi ■ k 



k 4 



d n l 



l 2 (k~l) 2 

Vi ■ I 



d n l- 



l 2 (k-lf 
(Vi ■ I) 2 



k 4 J l 2 (k-l) 2 
Expanding around Vi ■ k' — this gives 



(50) 



ji 

1 3bgl 



xT 



if 2- 5+2e i7r- 2 +^(l - /3 2 ) 



r i 



r(2-e) 



(10 - 4e) 1 



d n k 



3-eJ J (Vi -fc) 2 (fc 2 ) 1+ f (51) 
Using Eq. ([TT|) the UV poles of the integral are then 
a 2 19 



L 3bgl 



ir 2 48 





58 TS 




{---- 


VKa 

57 " 


} 


\ e 2 e 





(52) 



The counterterm is I^ gl = (a 2 /7r 2 )(19/48)[2/, 
Kp/e}. Then 

a 2 19 f 1 58 
^48 1 T 2 ~ 57e 



'3bg2 



(53) 



Last we calculate the ghost-loop self-energy dia- 
gram in Fig. 3(b). After several manipulations and 
expanding around Vi ■ k' — this gives 

I'^ h = i 2 - 6+2£ 7r- 2+ *z(l - (3 2 )T (-1 + |) 



<3bgh - ^2' 



1 



d n k 



T(A-e)J (vi ■ k) 2 (k 2 ) 1+ i 
Using Eq. ([77|) the UV poles are then given by 



i 36 9 h ^2 48 1 t 2 



1 1 

e 2 ~ ~e 

c.t. _ (-^,2 /^2 



4 

3+^ 



.(54) 



(55) 



with counterterm I^ gh = (a 2 /ir 2 )(l/A8)[2/e 2 + Kp/e 
Then 



a 2 1 



'3bgh 



7T 2 48 e 



1 



4 

3e 



(56) 



The sum of the gluon and ghost loops, Eqs. ([5 
and ([56"]) . denoted by hbg, is then 



, a? 5 

i3bg - — — 



7T 2 12 



1 



31 

We 



(57) 



G 
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4.3. I 



3c 



Finally we calculate the diagram in Fig. 3(c). Using 
Eqs. ([69]) and ([70]) and adding the one- loop countert- 
erm, we find 



3c 



a2(l + /3 2 )(-l), fl-0 



2/9 



In 



1 + /3 



(58) 



5. Two-loop soft anomalous dimension 

We now combine the kinematic results from sections 
3 and 4 with color and symmetry factors. The con- 
tribution of the diagrams in Figs. (2) and (3) to the 
two-loop soft anomalous dimension is 



Cp [l2a 

f Cf Ca 
1 



- -t- 
-\h b 



2 I 2 d + 2 l2e - 
+ hf~ heg 



'3cJ 



:C F [I; 



2cq 



l 3bq\ 



2e 2 V s J 4e 2 s 2e 



(59) 



On the right-hand side of Eq. ([59]) in addition to 

(2) 

Tg, which appears in the coefficient of the 1/e pole, 
there also appear terms from the exponentiation of the 
one-loop result and the running of the coupling, with 
[3q — (11/3)CU — 2nj/3, Ca = N c , which account for 
all the double poles of the graphs. From Eq. ([5"9"]) we 
solve for the two-loop soft anomalous dimension: 



r 



(2) 



K 



2 6 

-/3 2 ) 2 



+ CfCa 



8/3 2 



Cs + C2 In 



2 



1 



In" 



1-/3 



In 



1-/3 



1 - 

/3 2 ) 



4/3 



Li 2 



-Cain 

/3 



1 + /3 

(1-/3) 2 
(1 + /5) 2 
1-/3 



1 + /3 



5" 



Li, 



In" 



1 + /3 
, / 1 - /8 



1 + /3 

(1-/3) 2 
(1 + /3) 2 

1-/3 



1 + 
- Li 2 



2 In 



1-/3 
1 + 



hi 



1 + /3 
(1 + /3) 2 



4/3 



where 



(60) 



(61) 



A slightly different but equivalent expression is given 
inRef. [|. 

In terms of the cusp angle 7 we find 



r 



(2) 



K 



(i) 



CfCa 



^1 
2 



7_ 
2 



1 



■ coth 7 





50 








40 








30 






y' 










20 






>* 








■/ 


10 
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f 2) 

Figure 5: The two-loop soft anomalous dimension Ig . 



C 3 - C27 - y ~ 7L12 (e^) - Li 3 (e" 27 ) 



- coth 7 



2 7 
C2 + C27 + 7 + y 



+ 2 7 In (1 - e" 3 ?) - Li 2 (e"^)] } . (62) 



This expression is consistent with the form of the 
two-loop cusp anomalous dimension given in Ref. Q 
(which involves a few uncalculated integrals), and it 
is also consistent with the two-loop heavy-quark form 
factor of Ref. [6j . 

Using the above results we can find the soft- 
gluon logarithms in the cross section for e + e~ — > tt, 
which are of the form ln" _1 (/3 2 )//3 2 at n-th order in 
a s . The first-order soft-gluon corrections are c' 1 ' = 
a B (a g /ir) 2 Tg/0 2 with a B the Born cross section. 
The second-order soft-gluon corrections are 



4(r«) 2 



/3 r 



(i) 



2TiT 



(i) 



2T 



(2) 



ln(/3 2 ) 

P 2 
1 



-7T2 (63) 



with T\ the NLO virtual corrections. 

In Figure 5 we plot Tg ' versus p. The insets show 
the small and large regions in detail. Note that both 



Vg and To 1 ' vanish in the threshold limit, = 0, and 
they both diverge in the massless limit, (3 = 1. 
The small- /3 expansion of Eq. ([50]) is 



.(2) 



.(2) 
S exp 



27 



/3 2 [C f Ca(18( 2 



47) + 5n f C F } 
0(/3 4 ). (64) 



(2) 

F s is an even function of and hence only even pow- 
ers of 13 appear in the expansion. As shown in [3| the 
expansion provides a good approximation to the com- 
plete result for small (3. 
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(2) 

Next we study the large-/? behavior of Fg , 



([50j). As /? 



r 



(2) 



A 



(1) 



CfCa 



(l-Ca) 



Eq. 



(65) 



In (2u 
diverges at (3 = 1, T 



with = c F 



1 



Since T 



(i) 



(2) ~ (K/2)T ( g ) 



at that limit. 

2) R-t^(I) 



This is consistent with the massless case T s = ^-T s 
with ri 1 "* = Cp ln(Uj 



/) As is clear from Eq. 

60| the massive case is much more complicated than 
the simple massless relation, and (K/2)T^ is just the 

(2) 

first of many terms in the expression for T s . Figure 6 

shows that numerically the ratio (K/2)T^ /T^ goes 
to 1 as (3 — > 1, the massless limit, as expected from the 
above discussion, but it is significantly different from 
1 at other (3 and takes the value 1.144 at the other 
end of the range, (3 = 0. 

In the mixed massive-massless case, with Vi the 
heavy quark and Vj the massless quark, we find 



r 



(2) 



K 
~2 



with 



lf> = C F 



r (l) 

1 s 


+ CfCa 


-GO 

4 


H 


y/2 Vi ■ Vj \ 


1 




v V^F ) 


~ 2 



(66) 



(67) 
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Figure 6: Approximations to I5 . 



next-to-next-to-leading logarithm (NNLL) resumma- 
tion for heavy quark hadroproduction and has been 
used in very recent calculations [3, EE Ell E3 • 
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Given that the small-/? expansion gives very good 
(2) 

approximations to T s at smaller p while the expres- 
sion (K/2) is the large-/3 limit, we can derive an 

approximation to T ( s , Eq. (f6"0")) . valid for all (3 using 
the following approximate formula: 



.(2) 
S approx 



F 



(2) 

S exp 



2 s 



(i) 



(i) 



CfCa 1 



A\ 



^ 2 



-(i) 

S cxp 

/3 2 + (/3 4 ) . (68) 



As seen from Fig. 6, adding just the 1 terms 
to (A/2) Tg provides an excellent approximation to 
the exact result for for all (3. If we keep additional 
terms through order (3 12 then the approximation is 
extremely good and does not differ by more than one 
per mille anywhere in the f3 range. 

Finally we discuss extensions to single top [l(| EH 
and top pair [l(| EH production at hadron colliders. 
As first shown in Ref. [H, since the two- loop soft 
anomalous dimensions for single top and top pair pro- 
cesses involve eikonal graphs with both massless and 
massive eikonal lines, it is clear that the simple mass- 
(A/2) ri 1} will not hold. This 



less relation r^ 2 -* 



observation was also made later in Refs. [H, [Til ]. 

(2) 

The two-loop soft anomalous dimension T y s ' derived 
in [H and in this paper is an essential ingredient for 



Appendix 

We list results for the UV poles of several of the 
many integrals needed in the calculation of the two- 
loop soft anomalous dimension. 



d n k 



k 2 Vi ■ k Vj ■ k 



l(-l)- 1 -i 7T 2 -i2 



3+3 f 



(3 \e \l + f3 
+ l(21n2 - lnvr - 7E - ztt) In 

+ iln 2 (l+/3)-lln 2 (l-/3) 

- \ hh (Hr ) + \ hh C-2 1 



0(e) 
(69) 



where 2-F1 is the Gauss hypergeometric function. 



d n k 



k 2 (vi ■ k) 2 



3+3-s 
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x(l-/? 2 )- 1- *r(l 



I) 



T^ 2 



1 

~ 2 
In 2 



ln(l-/3 2 ) 



• ln7r 



IE 

2 



J7T 

y 



•O(c).(70) 



cPfc 



(fc 2 ) 1+ * (Uj • fc) 2 



x(l-/3 2 )- 1 - £ r(l + e) 



^- e 2 2+3e 7T 2 - 
1 



r i 



(77) 



d n k 



i (-iy 



l2e 



X7r 2 -*r(l + e) 



r(i + i) 



{l-!3)-\F l ^- e ,l + e;l-e; 
- (l+f3)- e 2 Fi (-e,l + e;l-e 



2 

1 + /3 



d n k 



k 2 {v t ■ k) 1+e Vj ■ k e(l + e) 



x(i-/? 2 r 



- r M 



3e 



2 2+ ~(-l) 



1 



-l— as 



2 y r(i + e) 



3e 3e „ 2/3 -2/3. 

x^[i + £ ;i + y ,i + Y ;2 + e ; T -^, T -^] (72) 

where F\ is the Appell hypergeometric function. 



d n k 2 



i (-1)"^ 



k\ [v. l -{k 1 + k 2 )f e (1 + e) 

x (i- ( 3 2 )- 1+ Hv i -k 1 )-<r(i + ^ 

3 + e 



xr(i--)r 



d n k 



(73) 



k 2 (vi • fc) 2 + e e(l + e) 



x(i-/3^)- i --r i + - 



3c 



1 



k\ Vi ■ ki Vi ■ (ki + k 2 ) 



2 J r(l + e) 



(-1)2 2 2 ~4tt ; 



(74) 



x(« 4 -fc 2 )- e (i-/3 2 )- 1+ *r(i + |) 
xrfi-l)^). (75, 



d n k 2 



k 2 Vi ■ k 2 [vi ■ (hi + k 2 )] 2 1- 



(-1) 



i-H 



x2 3 -^ tt 2 -4 («j • hy 1 -" (1 - /3 2 )" 



xr(i--j r(i + e ). 



(76) 
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